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Abstract 


We evaluate in detail the properties of nuclear gamma-ray emission in the 
4 to 8 MeV range. This emission consists of broad (AE/E ~ 20%) and narrow 
(AE/E ~ 1%) lines resulting from nuclear reactions of energetic H, He, C 
and 0 nuclei with ambient matter. From the comparison of the calculations 
with observations of the 1972, August 4 flare we conclude: 1) that essen- 
tially all the observed radiation in the 4 to 8 MeV region is due to the 
superposition of broad and narrow lines of nuclear origin with almost 
no contribution from other mechanisms; 2) that the accelerated particles in 
the energy region from about 10 to 100 MeV/amu have a relatively flat 
energy spectrum (for a power law, E , s^2.5); 3) that the calculated 
gamma-ray spectrum obtained from an isotropic distribution of accelerated 
particles fits the observed spectrum better than the spectrum derived from 
an anisotropic distribution for which the particles' velocity vectors point 
towards the photosphere; and 4) that it is possible to set a stringent 
upper limit on the ratio of relativistic electrons to protons in flares, 
consistent with the small, but finite, electron-to-proton ratio in galactic 


cosmic rays. 
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I . Introduc tion 

The solar gamma-ray spectrum from about 4 to 8’ MeV consists of 
strong nuclear lines superimposed on continuum emission which could 
arise both from nuclear interactions and bremsstrahlung of energetic 
electrons. These nuclear reactions are almost exclusively due to 
interactions between hydrogen, helium, carbon, nitrogen and oxygen 
nuclei. Reactions due to energetic protons or alpha particles produce 
narrow lines with full widths at half maximum (FWHM) less than 150 keV; 
reactions caused by energetic C, N and 0 nuclei give rise to broad 
lines which in most cases cannot be distinguished from the continuum. 

Line emission at 0.51, 2.2, 4.4 and 6.2 MeV and continuum emission 
extending up to about 7 MeV associated with the large solar flare of 
1972, August 4 were reported by Chupp et al. (1973, 1975). These lines 
are due to positron annihilation, neutron capture by protons ? and 
deexcitation of excited states in C and 0 (e.g. Ramaty, Kozlovsky and 
Lingenfelter, 1975). Electron bremmstrahlung could account for & significant 
fraction of the continuum, at least up to energies of several MeV (Suri 
et al. 1975; Bai and Ramaty 1976). 

The production of gamma rays in solar flares has previously been 
studied by several authors (e.g. Lingenfelter and Ramaty 1967, Ramaty, 
Kozlovsky and Lingenfelter 1975, Crannell et al. 1976). Even though 
these papers provided detailed studies of the 0.51 and 2.2 MeV lines, 
they limited their treatments to only the gross features of gamma-ray 
production in the 4 to 8 MeV range. In particular, the shapes of the 
gamma-ray lines and the origin of the continuum in this energy band have not 
been investigated in detail. Ramaty and Crannell (1976) have provided 
a preliminary study of the shape of the 6.13-MeV line of ^0 and showed that 
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the mean energy of this line depends on the angular distribution of the 
primary energetic charged particles. They have limited their treatment, 
however, to the narrow line component only, whereas, as we shall show in the 
present paper, very important energy shifts do take place also in the broad 
components . 

The purpose of the present paper is to provide a detailed treatment of 
nuclear gamma-ray production in the 4 to 8 MeV region. This energy band 
is of particular interest because it contains essentially all the lines of 
the abundant nuclear radiators, C and 0, and conversely, it is not signifi- 
cantly populated by photons from other nuclei. Furthermore, as we shall 
show, the entire observed radiation in the 4 to 8 MeV range could be due 
to the superposition of broad and narrow lines from C and 0 with essentially 
no contribution from other radiation processes. This result allows us to 
carry out, for the first time, a detailed comparison between theoretical 
and astronomical gamma-ray line spectra. 

In section II we examine in detail the nuclear cross sections. We find 
that the extant literature contains sufficient data to present a self- 
consistent picture for the reactions leading to gamma-ray production in 
the 4 to 8 MeV energy band. However, several key questions remain to be 
resolved. In particular, future experimental work, using existing laboratory 
accelerators and high resolution gamma-ray detectors should be able to resolve 
all the major lines in this energy band, and to provide information on the 
cross sections of several lines which have not been measured in past experi- 
ments using low resolution instruments. In Section III we use the nuclear 
cross sections to calculate the production of gamma rays, and we investigate 
the dependence of the shape of the gamma-ray lines on the energy spectrum and 
angular distribution of the primary accelerated charged particles. We use 
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detailed relativistic kinematics, and we take into account, whenever 
available, data on the angular distribution of the recoiling secondary 
particles . In Section IV we compare the results of our calculations 
with the observational data; in particular we present a direct comparison 
between the observed pulse-height spectrum and the calculated gamma-ray 
spectra folded through the response function of the University of New 
Hampshire gamma-ray detector flown on 0S0-7. We apply the results of 
our calculations to solar flares, because until now no gamma-ray lines 
were convincingly observed from other astronomical sources. However, 
our results and techniques are relevant to future gamma-ray line observa- 
tions from outside the solar system as well. We summarize our results 
in Section V. 
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II. Nuclear Cross Sections 

The principal gamma-ray lines in the energy range from 4 to 8 MeV, 
their emission mechanisms and production modes are summarized in Table 1 0 

Line emission at 4.44 MeV is produced by the deexcitation of the 
12 

first excited state of C at 4.439 MeV and of the second excited state of 

H-B at 4.444 MeV. Because of the recoil of the nuclei emitting gamma rays, 

the energies of the photons are 4.438 and 4.443 MeV respectively 0 The 

i 2 12 A A on 

cross section for the reaction C(p,p' ) C* MeV has been measured 

by Reich et al. (1956) between 5 and 5.7 MeV, by Barnard et al. (1966) 
between 6 and 11 0 5 MeV, by Conzett (1957) between 10 and 12 MeV, by 
Daehnick and Sherr (1964) between 14 and 19 MeV, by Dickens et al. (1963) 
between 18 and 30 MeV, by Stovall and Hintz (1964) at 40 MeV, by Fannon 
et al. (1967) at 49.5 MeV, by Horowitz and Bell (1970) at 100 MeV, by 

Emmerson et al D (1966) at 145 MeV, and by Tyren and Maris (1957) at 

185 MeV. This cross section is shown by the solid line in Figure 1. We 
have used the average cross section at the resonances near the maximum 
at 10 MeV. However, we have not averaged the cross section at the resonance 
at 5.35 MeV since it may have some effect on the calculations for very steep 
particle spectra. The data points represented by the closed circles and 
square are based on measurements of gamma rays at 4.44 MeV resulting from 
the bombardment of C with protons (Zobel et al. 1968, Alard et al. 1974). 

The closed circle points were deduced by using the differential cross 
section at 135° given by Zobel et al. (1968) and their data on the angular 
distributions of gamma rays produced in proton bombardment of Co The data point 
of Alard et al. (1974) is already given in integral form. As can be seen, 
the gamma-ray data is consistent with the (p,p' ) data represented by the 
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solid line up to ~ 20 MeV. The discrepancy at higher energies is most likely 

due to the excitation of the 4.444 MeV level of B by the reaction 

12 C(p,2p) 1X B which has an energy threshold of 22 0 1 MeV„ This discrepancy 

12 

cannot be due to the excitation of higher levels of C which cascade to 

12 

the 4.439 level because all the levels of C above this level decay almost 

exclusively by particle emission (Ajzenberg-Selove and Lauritsen 1968) 0 

Since the kinematical widths of both 4.438 MeV and 4.443 MeV lines are 

~ 50 keV, it is impossible to distinguish between these two lines. Therefore, 

in our calculations for the production of the 4.44 MeV line above 

20 MeV we use the cross sections given by the dashed line in Figure 2. 

Up to about 100 MeV these cross sections are higher by approximately 

25% than those used by Ramaty et al. (1975) who took into account only the 

(p,p # ) data, but they are lower than the cross sections given by Meneguzzi 

and Reeves (1975), because these authors assumed complete isotropy 

for the gamma-ray data of Zobel et al. (1968). 

12 12 *4 0 44 MeV 

The cross sections for the reaction C(or t </ ) C are given 

by the dashed-dotted line in Figure 1. The measurements of these cross 
sections were summarized in Ramaty et al 0 (1975). The open circle represents 
a measurement by Zobel et al. (1968) for 4.44— MeV gamma ray production 
In ct-C interactions,, As in the case of proton induced interactions, we 
interpret the difference between the gamma-ray data and the (&,</ ) data 
as due to the contribution of ^Bo In our calculations, we use the 
dashed-dotted line up to 5 MeV/ ami and the dashed-crossed line at 

higher energies. 

12 11 

The 4.439 MeV level in C and the 4.444-MeV level in B can also 
be populated by spallation reactions on ^0, as indicated in Table 1. 
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The cross sections for these reactions were taken from Zobel et al. (1968) 

and Alard et al 0 (1974), and are shown by the solid lines in Figure 2. For 

the (a, 2a) reaction we have only one experimental point at 13 MeV/nucleon. 

Our estimate of the cross section curve for this reaction la based on the 

expected behavior of the cross section near the threshold and an arbitrary 

extrapolation to higher energies where the contribution of this reaction 

is negligible in comparison with the proton induced reaction. Spallation 
14 

reactions on N can also lead to 4.44-MeV gamma-ray production. Using the 

measurements of Macleod and Reid (1966) at 13 MeV and those of Clegg et al. 

(1961) at 120 MeV, we estimate that the 4.44-MeV gamma-ray cross section 

from is larger than that from ^0 by about a factor of 2. With 

14 

solar abundances, N/0 ~ 0,18, the 4.44-MeV yield from N is about 30% 

16 

of the yield from 0 o 

Nuclear excitation of ^0 can produce several gamma-ray lines. The 
16 

first excited state of 0 at 6.05 MeV decays by e+ ■>- e~ pair emission. 

The second, third and fourth excited states at 6.131, 6.919 and 7.119 MeV 

decay almost exclusively to the ground state producing gamma rays at 

6.129, 6.917 and 7.117 MeV, respectively. The fifth excited state at 

8.872 MeV decays 75% of the time to 6.131 MeV level and hence produces 

lines at 2.741 MeV and 6.129 MeV. Other transitions between the states 
16 

of 0 are also possible (De Meijer, Plendl and Holub 1974). In 
particular, excited levels at 10.94 and 11.07-MeV cascade to the ground 
state via the excited states at 6.131, 6.919 and 7.119 MeV with branching 
ratios given by Aj zenberg-Selove and Lauritsen (1968). 

1 C 

Gamma-ray lines from '*' u '0 can also be produced by spallation reactions. 
The reaction ^0(p,2p)^N and ^ 0 (p,pn )^0 can populate excited states in 
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and ^0 o The strongest lines from are at 6.322, 5.270 and 5.298 MeV 

and from ^“*0 at 5.180, 5.241 and 6.176 MeV. 

The cross section for producing gamma rays at 6.129 MeV from proton 
16 

interaction with 0 is shown in Figure 3 by the solid line. This 
cross section is sum of the cross sections of the reaction ^OCpjp' ) 

16q* 6.131 Q f cross sec tion for the reaction ‘^0(p,p / ) ^0 

and about 50% of the cross section of the reaction ^^0(p,p / ) ^^0 These 

cross sections were measured by Dangle et al. (1964) between 7 and 10.5 MeV, 

by Kobayashi (I960) between 11 and 15 MeV, by Daehnick (1964) between 15 

and 19 MeV, by Hornyak and Sherr (1955) at 19 MeV, by Crawley and Garvey (1967) 

at 17.5 MeV, by Austin et al. (1971) between 17 and 45 MeV, and by 

Sundberg and Tibell (1969) at 185 MeV. Since the contributions of 

cascades from higher levels to the 6.131-MeV level was neglected by 

Lingenfelter and Ramaty (1967) and Ramaty, Kozlovsky and Lingenfelter 

(1975), the cross sections for 6.129-MeV line emission used by these 

authors were somewhat lower (~ 30-50% between 10 and 50 MeV) than those 

used in the present calculation,, 

There are no measurements of 6,129-MeV line emission from spallation 
20 

reactions on Ne. If we assume that the cross section for the reaction 

20 .16 *6ol29 , _ , . 16^, . 12, *4 o 439 

Ne(p,pa) 0 is similar to that of the reaction 0(p,pa) C 

20 

then, with solar abundances, the contribution of Ne to 6.129-MeV emission 

X6 

is less than about 10% of the contribution of 0 o 

The data points shown in Figure 3 are based on the measurements of 
Zobel et al. (1968), Foley et al. (1962), and Alard et al. (1974). In 
these measurements the cross sections were determined by measuring gamma 
rays at about 6 Q 2 MeV. Because of the poor energy resolution of the 
detectors used, these gamma rays include the lines at 6.129 MeV from 



^0, at 6.176 from ^0 and 6.322 MeV from ^N, The difference between 

the dashed and solid lines in Figure 3 gives the sum of the cross sections 

for the production of the 6.176 and 6.322 MeV lines. This cross section 

is shown by the dashed-dotted line in Figure 2. There is only one 

16 

direct high-resolution measurement of gamma rays from 0 which can 
resolve the 6.129, 6.176 and 6.322-MeV lines (Goryachev et al. 1973). 
According to this measurement the 6.322-MeV line is stronger by about 
a factor of 2 than the line at 6.176 MeVo Assuming that this ratio holds 
also at lower energies, we take the cross sections for the production of 
the 6.322-MeV and the 6.176-MeV lines as 2/3 and 1/3, respectively, of the 
dashed-dotted curve in Figure 2. The sum of the cross sections for the reac- 

i r *i c g I z: I c ^ 

tions 0(p,pn) 0 and 0(p,2p) N leading to ^ 5.2-MeV photons 

is shown by the dashed line in Figure 2 based on data by Zobel et al D (1968) 

Most of the experimental measurements of the cross sections for the 
excitation of the 6.919 and 7.119 MeV levels in 0 by (p,p' ) reactions do 
not resolve these two levels. Therefore, in Figure 4 the solid line 
gives a fit to the sum of the cross sections for the excitation of these 
levels based on measurements of Zobel et al 0 (1968). Measurements of 
(p,p') reactions to these levels by Kobayashi (1960) at 5.6 MeV and 
Crawley and Garvey (1967) at 17.5 MeV coincide with this line. At higher 
energies the measurements of (p,p') reactions by Hornyak and Sherr (1955) 
and Sundberg and Tibell (1969) give lower values. The difference should be 
attributed to contributions of spallation lines such as the 7.3 MeV from 
^~*N, which show up at high energies. 

The 2*741 MeV line is produced by transitions between the 8.87 MeV level 
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to the 6.13 MeV level in ^0 (branching ratio 76%). The dashed line in Figure 4 

is based on the measurements of Crawley and Garvey (1967) at 17 <,3 MeV, and 

of Austin et al. (1971) between 17 and 45 MeV. The data points are from 

Zobel et al. (1968). These data are in agreement with the (p,p' ) measurements. 
16 

The a + 0 -* Y , no ._ cross section is shown by the dashed-dotted line 

b.y Z+/ , 

in Figure 4. This line is based on measurements of Corelli et al (1959) at 
4.5 MeV/amu, of Harvey et al. (1966) at 10 MeV/amu and of Zobel et al. at 
13 MeV/amu. 

In the interval 4 to 8 MeV the .n-eractions of protons with C produce a 

cluster of lines around 6.5 MeV due to deexcitations of excited states in 
11 11 

C and B. These are populated by (p,pn) and (p,2p) reactions which are 

known to be important processes leading to spallation. The cross section 

for this cluster, as given by Zobel et al. (1968) and Clegg et al. (1961), 

is on the average about 1/4 of the cross section for the production of the 

16 

~ 5.2-MeV feature from 0 by the same type of reactions. We use this 
ratio in our subsequent calculations 
III. Gamma-Ray Production 

.1. Total Narrow Line Production Rates 

The instantaneous reaction rates between energetic protons or 
alpha particles and the ambient medium for the processes listed in Table 1 
can be written as 

q, ± (p) = n ± J* dE CTi (E)cg N p (E)dE (sec" 1 ) (1) 

and 

q.(or) = n. ? dE a . (E)cp N (E)dE, (sec -1 ), 

n X X X O' 

where n. and <j . are the target number density and cross section for 


( 2 ) 
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the particular reaction considered, E and c[3 are energy per amu and velocity 
of the projectiles, and N^(E)dE and N^(E)dE are tne total number of protons 
and alpha particles in the interaction region with energies per amu in dE 
around E. The instantaneous rates given by equations (1) and (2) are valid 
for both the thin and thick target models (e.g. Ramaty et al. 1975). The 
relationships between and N^ and the source functions of these particles 
are different, however, in the two models. 

We have evaluated equations (1) and (2) for the cross sections given 

' -3 

in Section II, ambient abundances normalized to 1 H atom cm by the ratios 
H:He:C:N:0 = loSlO^:3. 7:1. 2:6.8 (Cameron 1973, except for the He/H ratio) 
and for proton and alpha particle spectra of the form 


N (E) = ION (E) = (s-l)30 S_1 
p a 


X' 


fE S E E 

c c 

E -s E £ E 

^ c 


(3) 


These spectra are normalized such that for E^ ^ 30 MeV/amu and s > 1, the 
integral of N^(E) above 30 MeV is unity. We use power-law spectra with low- 
energy cutoffs, because, as we shall see below, significant changes in the 
gamma-ray spectra are caused by variations in the value of E c< 

The resultant production rates of narrow gamma-ray lines at the 
indicated energies are shown in Figure 5 as functions of the spectral index 
s and the cutoff energy E c « The production rates of these narrow lines 
are due to the interaction of energetic protons and alpha particles with 
ambient C, N and 0. The inverse reactions induced by energetic C, N and 0 
nuclei leading to broad lines are not included in the curves in Figure 5. 

The 4.44 curve of Figure 5 represents the sum of the 4.438 and 4.444~MeV 
lines listed in Table 1; the 5.2 curve is the sum of the 5.180, 5.298, 

5.270 and 5.298 lines; the 6.13 curve represents the 6.129-MeV line; the 
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6.2 curve is the sum of the 6.176 and 6.322-MeV lines, but as discussed in 
Section IX, the 6.322 line appears to be stronger than the line at 6.176 MeV; 
the 6.5 curve is the sum of the 6.337, 6.478, 6.741 and 6.791-MeV lines; 
and the 7 line is the sum of the 6.917 and 7.117-MeV lines. 

As can be seen from Figure 5, the strongest line in the 4 to 8 MeV 
range is at 4.44 MeV. The ratio between its intensity and that of the 
next strongest line, at 6.13 MeV, is about 2, independent of s or E^. The 
line at 6.32 MeV and the feature at ~ 5.2 MeV, due to spallation reactions, 
are relatively strong for flat particle spectra, but become quite insignificant 
for accelerated particles with steep energy spectra. 

The curves labelled and q^ represent the sum of the production 
rates of the narrow lines listed in Table 1, due to proton and alpha 
particle interactions, q^ = £ q^(p) and q^ = Sq^or) . As can be seen, 

9p + 9cy equals to the sum of the production rates of the various lines 
shown in this figure. 

2. Gamma-Ray Spectra 

We proceed now to calculate the spectrum of gamma rays resulting 
from the above interactions by taking into account in detail the kinematics 
of the interactions. We first discuss the narrow gamma-ray lines resulting 
from the direct excitation reactions. 

Consider the reaction between a proton or alpha particle of mass m^ 

and energy per amu E, and a target of mass m^ which is stationary in the 

frame of the Sun. This reaction leads to an excited nucleus of mass 

m„ = m„ + e , where e is the energy of the excited level. The Lorentz 
3 2 x x 

factor (total energy per rest mass energy) of the excited nucleus in 
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the center-of-mass (c.m.) frame is given by 


2 2 2 

Y* = (E + m_ - m 1 )/(2m_E ), 

> c rm ^ 13 cm 


1 s 


cm 


( 4 ) 


where 


(5) 


2 i 

E = [(m.+nu) + 2m .. itu E/ml 2 

cm 1 2 12 p 

is the total energy in the c.m. frame and m^ is the proton mass. In 

terms of this Lorentz factor and the directional cosine of the excited 

nucleus in the c.m,. frame, u, , the Lorentz factor and directional cosine 


of the excited nucleus in the frame of the Sun are given by 


= y c Y s 


( 6 ) 


*- Y s y c 2 -1 + y c v / y* 2 -1 


(7) 


where 


( 8 ) 


(9) 


y = (E 2 + m 2 -m 2 ) /(2m.E ) 

Y c cm l 2 1 cm 

is the Lorentz factor of the c.m. 

We denote by u ' the cosine of the angle between the gamma ray 

Y 

and the velocity vector of the excited nucleus in the rest frame of 
this nucleus o In the frame of the Sun this directional cosine is 

= (P s + t V' )/(1 + 

where B = Vv 2 -i/v . Taking into account the Doppler effect and the 
s , s 1 s 

recoil of the nucleus during photon emission, the gamma ray energy 
is given by 

£ = A£[1-A£/(2Mc 2 )]/[y (1-P u ) ] , 

Y soy 

where A£ is the difference between the energies of the initial and final 

states of the transition and M is the mass of the gamma-ray emitting nucleus 

Equation 10 is valid if the gamma ray is emitted while the nucleus is still 


( 10 ) 




in flight. Only in solids will some nuclei stop before emitting photons, ine 
astrophysical consequences of this possibility were discussed by Lingenf elter and 
Ramaty (1976) . Finally, the cosine of the angle between the gamma ray 
and the velocity vector of the projectile is 






where d> -d> is the difference between the azimuths of the projectile and 
s y 

gamma ray. 

We use a Monte-Carlo simulation to calculate the energy distribution 
of the gamma rays. The distribution of E is determined by the total cross 
section as a function of energy and the energy spectrum of the primary 
particles . 

The distribution of p, g * is determined by the differential cross 

section in the c 0 m. frame. At low energies 10 MeV) , where resonances 

occur, there are rapid variations of the angular shapes, but on the 

average the distributions are close to isotropic. At higher energies, 

the differential cross sections are more and more peaked in the forward 

direction, and hence the excited nuclei are emitted preferentially backward 

with respect to the direction of the incoming protons. In order to 

account for these effects, we evaluate the differential cross sections 

in 4 different regions of energy. Below 12 MeV we assume an isotropic 

distribution; between 12 and 20 MeV we use the data of Peelle (1957) at 

15 MeV; between 20 and 60 MeV we adopt the measurements of Dickens et al. 

(1963) at 28 MeV; and above 60 MeV we use the measurements of Strauch and 

Titus (1956) at 96 MeVo All these data are for 12 C(p,p' ) 12 C reactions. 

16 16 *6 13 

By examining the corresponding data for the reaction 0(p,p) 0 
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(Kobayashi 1960, Daehnick 1964, and Austin et al. 1971) we conclude 
that, on the average, the shape of the differential cross section of this 
reaction is almost the same as that of the ^^C(p,p / )^C reaction. 

By using a random number generator to evaluate E and p g *, we can 
now calculate, from equations (4), (5), (6) and (8), the energy of the 
excited nucleus. The energy of the gamma ray then follows from equations 
(9) and (10) , where we assume an isotropic distribution of the gamma rays 
in the rest frame of the emitting nucleus. The cosine of the angle 
between the gamma ray and the projectile, p, can now be calculated by using 
equations (7) and (11) , where we assume a uniformly distributed random number 

for 4> -<J> . The knowledge of p is required only if the projectile velocities 

s y 

are anisotropic. For isotropic distributions all the gamma rays generated 
by the simulation contribute to the spectrum, regardless of the 
value of p. However, by imposing conditions on p the effects of anisotropic 
distributions can be taken into account. For example if the projectile 
velocities are uniformly distributed in the hemisphere away from the 
observer, only gamma rays with -l^p^O contribute to the spectrum, while if 
the projectiles are confined to the opposite hemisphere, (Kp^l. Hereafter 
these anisotropic distributions will be referred to as the down and up 
cases, respectively. 

In calculating the shape of the gamma-ray lines from the spallation processes 
of Table 1, we use the same techniques as for the direct excitation processes, 
except for the following modifications. The Lorentz factor y s * is no longer 
a unique function of the incident energy E, but can have any value up to 


a maximum 


V * 


* = (E 2 -f m 2 -m 2 )/(2 m-E ), 
i cm 3 4 3 cm 


( 12 ) 
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where m^ is now the mass, including the excitation energy, of the excited 

nucleus and m^ is the sum of the masses of all the other particles that 

emerge from the reaction. Since there is no data on the distributions 

of y g * and jj, g * for these reactions, we have made certain assumptions which 

lead to a gamma-ray spectrum consistent with the observations of Goryachev 

15 

et al. (1974) who have resolved the 6.32-MeV line of N in the bombardment 

of ^0 with 1 GeV protons. The relatively narrow observed width of this 

line in the 1 GeV experiment cannot be due to gamma-ray emission from 

nuclei stopped by the target material, because the lifetime of the 
15 -14 

6.32-MeV level in N is very short (<4x10 sec, De Meijer, Plendl and 

Holub 1974). The observations of Goryachev et al. (1974), therefore, imply 

that the, nuclei are produced with very low velocities in the laboratory 

frame. This behavior can be simulated by using the same angular distributions 

for the multibody spallation reactions as for the 2-body direct excitation 

processes and by imposing the additional constraints that 

v * = Y * fo r y * < v and y * = Y for y * ^ Y • Similar assumptions 

T s Y m T s T c T s 'c T s 'c 

were made by Ramaty and Lingenfelter (1969) and Ramaty and Kozlovsky (1974) 

3 

for He production in multibody processes resulting from proton-alpha particle 

3 

reactions in which very low energy He nuclei were produced. 

In addition to the narrow lines discussed above, broad gamma-ray lines 
are produced by the interactions of energetic C, N and G nuclei with 
ambient H and He. For the evaluation of the shapes of these lines we use 
the same formulae as for the narrow lines with some modifications 0 In the 
2-body processes m^ is replaced by m 2 in equation (4) , and p, g * is replaced 
-p, * in equations (6) and (7) ; in the 3-body processes the excited nucleus is 
assumed to move with the velocity of the center of mass, P c = Jy ^-1/y ^ . 
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Using the above kinematics and nuclear data, and assuming that the H:He:C:N:0 
ratios of the energetic particles are the same as those of the ambient medium, 
we have evaluated the gamma-ray spectrum due to both the narrow and broad 
components. In Figure 6 we show the results for isotropic particle 
distributions with energy spectra characterized by (s=4, E c =0) , (s=2, E c =0) 
and (s=2, E c =100 MeV/amu) 0 As can be seen, for steep particle spectra (s=4) 
the most prominent lines are at 4.44 and 6.13 MeV, while for flat spectra 
(s=2), the spallation lines at 6.32 and ~ 5.2 MeV also become quite noticeable. 

The full widths at half maximum (FWHM) of the 4.44 and 6.13 MeV lines 
are summarized in Table 2. For s=4 the FWHM' s are significantly larger than 
for the other cases because for steep particle spectra the lines are produced 
predominantly in alpha particle interactions in which the excited nuclei have 
larger recoil velocities than in proton interactions. The FWHM's, however, 
do not depend strongly on the energies of the primary particles. For s=2 
the width of the 4.44 MeV line, for example, is essentially independent of 
the cutoff energy E c » For all three values of E c shown in Table 2, this line 
is predominantly due to proton induced interactions. As the energy 
of the primary proton increases, the c.m. angular distributions of the 
excited nuclei become more and more peaked in the backward direction. 

As a result, these nuclei acquire only small recoil velocities even if 
the energies of the primary protons are large. A similar result would hold 
for the 6.129-MeV line, except that for E c = 100 MeV, this line is 
effectively broadened into a feature of FWHM =- 120 keV by the contribution 
of the 6.17-MeV spallation line. 
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Another noteworthy feature of the gamma-ray spectra from alpha particle 
induced interactions (s=4 in Figure 6) is the presence of broad features 
in addition to narrow lines. In particular, the feature at 4.4 MeV, due 
to the interactions of energetic C nuclei with ambient helium, has a FWHM 
of only 500 keV, and hence can be considered as an identifiable broad 
gamma-ray line. The absence of such a feature in gamma-ray spectra 
resulting from flat particle spectra is due to the higher average energies at 
which proton interactions take place, and due to the fact that in interactions 
with protons the energetic nuclei suffer less energy loss than in interactions 
with helium nuclei. 

In Figure 7 we show results for anisotropic particle distributions. As 
pointed out by Ramaty and Crannell (1976), if the charged particles are 
confined to a hemisphere either away (down) or toward (up) the observer 
there should be a Doppler shift of the mean energy of the lines. We also 
obtain such shifts; however, since we have taken into account the angular 
distributions of the excited nuclei in the c.m. frames of the reactions, 
we find that the shifts are smaller than those given by Ramaty and Crannell 
(1976) who have assumed isotropic angular distributions. For the 6.129 MeV 
line, we find in the down case red shifts of about 15 to 20 keV for 
accelerated particles with s=2 and a shift of about 30 keV for s=4. For the 
4.44 MeV line we obtain red shifts smaller by almost a factor of 2. We 
conclude that such shifts can only be detected with high resolution 
detectors and sufficient counting statistics. 

In addition to the shifts of the mean energies of the narrow lines, 
anisotropic particle distributions also cause detectable differences 
in the shapes of the broad lines which appear as underlying continuum 
in Figure 7. As we shall see in Section IV below, these differences can be 
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studied by comparing the calculated and observed pulse-height spectra 
from the solar flare of 1972, August 4. 

IV. Discussion and Comparisons with Data 

We proceed now to compare the results of our calculations with 
observational data from the flare of 1972, August 4. Gamma-ray lines 
were detected from this flare by Chupp et al. (1973, 1975), The detector was 
flown on 0S0-7 satellite and has been described in detail by Higbie et al. (1973). 

1. Comparisons of Pulse-Height Spectra 

Figures 8 and 9 show comparisons of the observed pulse-height 
spectrum in the energy range from 4 to 8 MeV with calculated pulse-height 
spectra o The solid lines represent the pulse-height spectrum from the solar flare 
obtained by subtracting the antisolar quadrant counting rates 
from those in the solar quadrant. The dashed lines represent the calculated 
pulse-height spectrum obtained by folding the gamma-ray spectra of Figures 
6 and 7 through the detector response. The folding technique is discussed 
in the Appendix. The absolute normalizations of the calculated spectra was 
determined by minimizing the values of x 2 given by 


= Ci (c i 0bs - c i calC)3/c i calc > 

where C.°^ S and C. ca ^ c are the observed and calculated count rates in 
r i 


(13) 


200 keV bins as shown in Figures 8 and 9. The resultant minimal X 2, s are 
also given in these figures. It is evident from Figure 8 that the calculated 
spectrum for s=4 gives a poor fit to the data. (The possibility that 
X 2 = 37.6 by chance is less than 0.001). Much better fits are provided 
by the spectra for s=2, E c =0 and s=2, E c =100 MeV. By considering 
Figures 6 and 8 we see that the poor fit for s =4 is caused by the excessive 
number of photons around 4.4 MeV. As discussed in Section III, these 
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photons are due to the broad feature at 4.4 MeV produced by (a,c/ ) reactions. 

As we shall see below, steep particle spectra can also be ruled out 

from considerations based on the ratio of the 4 to 8 MeV gamma ray intensity 

to the 2.2-MeV line intensity. 

In Figure 9 we compare the data with calculated spectra from particles 
with different angular distributions. We see that the isotropic and 
up cases give a better fit to the data than does the down case. The isotropic 
and up cases are consistent with models wherein the charged particles are 
accelerated in a second stage mechanism by upward moving shock waves (Wild, Smerd 

and Weiss 1963, Frost and Dennis 1971, Sturrock 1974, BaiandRamaty 1976) , whereas 
the down case would correspond to models in which the gamma rays are produced 
by accelerated particles moving downwards into the photosphere D Unfortunately 
the data is not sufficiently accurate to completely rule out models in 
which the particles are beamed downwards. 

A further examination of Figures 8 and 9 reveals that the best 
agreements between the observed and calculated spectra are obtained with 
particle spectra that are flat at low energies (E c =30 and 100 MeV/amu). For 
such spectra there is a large contirbution from spallation reactions which 
produce the strong line at 6.32 MeV. This line can account, in part, for 
the observed counts in the 6.2 to 6.4 MeV bin. We should caution, however, 
that at least some of the observed counts in this bin could be due to 
6 to 6o2 MeV photons which were shifted into the next higher bin due to 
the nonlinearity of the relationship between photon energy and pulse height, 
and perhaps due to a slight gain shift of the detector’s electronic system. 
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2. Total 4 to 8 HeV Photon Fluxes 

We proceed now to evaluate the implications of the total number 

of photons_in the 4 to 8 MeV region. As discussed above, these photons 

are produced both as narrow lines due to proton and alpha particle 

interactions, and as broad lines caused by the reactions of energetic 

C, N and 0 nuclei. For the compositions used above (H:He:C:N:0 = 

4 3 

10 :10 :3. 7:1. 2:6. 8 for both accelerated particles and ambient medium) 
the total 4 to 8 MeV photon production rate is given by 

Q 4-8 = 2( W* (14) 

where q and q are plotted in Figure 5. However, some of the solar abun- 

p O' 

dances, in particular that of He in both the energetic particles and the 
ambient medium, are not very well known. If we take these abundances as free 
parameters, then equation (14) is modified as follows 

Q, Q = q [ 1+20 . 4 (ck’/p) ] +10q [ 1+20. 4 (He/H) ] (a/p) , (15) 

H-O p Oi 

where a/p and He/H are the He to H ratios in the energetic particles 
and the ambient medium, respectively. It is also assumed in equation (15) 

that in the energetic particles the C to N to 0 ratios are the same 
as those of the ambient medium and that He/(C+N+0) = 42 (Webber et al. 1975) „ 
We have evaluated equation (15) for several values of He/H and a/p. 

The range 0 o 07 < He/H ^ 0.1 is consistent with recent measurements (Brown 
and Lockman 1975), while the range 0.02 ^ a/p ^ 0.1 takes into account 
most of the observed variation of this ratio for solar flares. In Figure 10 
we plot Q^_g as a function of s and E c for He/H = 0.1 and a/p = 0.05, and 
in Table 3 we show the variation of Q^_g with He/H and a/p for two extreme 
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particle spectra. From Table 3 we see that Q^_g can vary by about a factor 
of 2 for very flat spectra, and by about a factor of 5 for very steep spectra. 
The significance of these variations will become apparent below, 
when we compare the ratio of the 4 to 8 MeV photon flux with the 2.2-MeV 
flux. 

The 2.2-MeV line is due to neutron capture by protons. Ramaty 

et al. (1975) have presented detailed calculations of neutron and 

2.2-MeV gamma-ray production in flares. Using their compiled cross 

sections, we have evaluated the production rate of neutrons, Q^, for the 

accelerated particle spectra and compositions discussed above. The 

results are given in Figure 11 for a/p = 0.05 and He/H = 0.1. The 

variation of Q r with or/p and He/ll is summarized in Table 3. The variation 

of Q n with or/p and He/H is of about the same magnitude as that of . 

The ratio g/Q n » therefore is almost independent of these parameters. 

From the neutron production rates, Q^, given in Figure 11 we evaluate 

the 2.2-MeV gamma-ray production rate, ^ by using the neutron-to-gamma 

efficiencies, f, given by Ramaty et alo (1975, figure 9). The ratio 

8 /Q 2 2 is plotted in Figure 12 as a function of s and E c for or/p = 0.05 

and He/H = 0.1. From the analysis of the Appendix, the observed 4 to 8 MeV 

-2 -1 

photon flux for the 1972, August 4 flare was 0.19 +0.02 photons cm s , 

-2 -1 

while the observed 2.2-MeV flux was 0.28 _j<).02 photons cm s (Chupp et al. 
1975). These values lead to an observed ratio of 0.68 +0.1 which is 
shown by the shaded area in Figure 12. 

By comparing the observed and calculated ratios of the 4-8 MeV flux 
to the 2.2-MeV line flux, we can obtain information on the spectrum of 
the charged particles. Before considering this comparison, however, let us 
discuss the principal sources of uncertainty in the calculations leading 
to the results of Figure 12. 
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The error in the nuclear cross sections can be estimated from the paper 
of Zobel et al. (1968) who give errors of less than 20% for the most important 
lines. Since we consider the sum of the lines in the 4 to 8 MeV range 
and since the total cross sections for most of these lines are taken from 
Zobel et aL (1968), we feel that the upper limit on the error due to 
nuclear cross sections probably does not exceed about 20%. 

The uncertainties resulting from compositional variation can be 
assessed from Table 3. We see that whereas for a given spectrum, g and 
Q n can both vary by up to a factor of 5, the ratio g/Q^ does not vary 
by more than 20% when or/p and He/H are allowed to change over their full 
assumed range of variability. Compositional uncertainties therefore, 
also do not lead to large errors in g/C^ 2 * 

Another source of error in the evaluation of Q„ _ is the abundance 

/ o / 

3 

of He in the photosphere (Wang and Ramaty 1974) 0 Our calculations are 
3 -5 

based on He/H = 5x10 , a value close to that measured in the solar wind. 

3 

By decreasing He/H arbitrarily to zero, g/C^ 2 is i- ncrease d by only 

less than a factor of 2. On the other hand, it is unlikely that the 

3 

photospheric abundance of He is much higher than the solar wind value, 

3 

since the He abundance has now been measured in prominences 

(Hall 1975), and its value there is close to that detected in the 

3 

solar wind. We conclude that the photospheric He abundance is also not 
a major source of error in the calculation of Q. ./Q 0 0 . 

Potentially, the largest source of error in the determination of 
Q, -,/Q„ 0 is the assumed angular distribution of the neutrons which produce 
the 2.2-MeV line. This angular distribution should resemble that of 
the primary charged particles. Our calculations of 2 are based on an 
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isotropic angular distribution. If the neutrons move preferentially downward 
toward the photosphere, the conversion efficiency, f, is increased by about 
a factor of 2, and Q. 0 /Q „ decreases accordingly. If, however, the 

neutrons move preferentially upward, f can be decreased by a large factor 
and hence g/^ 2 is f- ncreasec * by the same factor. 

In summary, we find that the values of g/Q^ 2 cou bd possibly be 
smaller than those given in Figure 12, but by not more than a factor 
of about 2. They could, however, be significantly larger, owing mainly 
to the possibility that the neutrons are preferentially moving away from 
the photosphere. 

By considering now the results of Figure 12, we see that steep particle 
spectra are inconsistent with the data c The same conclusion has been 
reached earlier from the comparison of the observed and calculated pulse- 
height spectra. This conclusion is unmitigated by the above uncertainties 
because the lowering of Q. 0 /Q 0 0 by about a factor of 2 allows only 
slightly steeper spectra to become consistent with the data. The increase of 
Q , „/Q „ due to the possible anisotropy of the neutrons leads to flatter 

spectra than those deduced from the nominal conditions of Figure 12. In 
fact, because of this uncertainty in Q^_g/Q 2 2 * most va l ues of s with 

= 100 MeV/amu are probably consistent with the observed 4-8' MeV to 2.2 MeV 
line ratio, and, furthermore, the values of s deduced for other E ' s should 
be considered as upper limits only 0 

For particle spectra with E^ ^ 10 MeV/amu, the best value of s is 

about 1.9, whereas for E =30 MeV/amu, s =- 2.5. We cannot deduce a value 

c 

for s for E c = 100 MeV/amu, but, as we have just discussed, this case 
can also be consistent with the data provided that the neutron-to-gamma 
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conversion efficiency is smaller by about a factor of 2 to 5 than the 
value given by Ramaty et al. (1975) for isotropic particle distributions. 

For particle spectra which remain unbroken power laws down to energies 
less than the gamma-ray production threshold energy (E = 0) , we can set 
a stringent upper limit on the value of s. If we allow a maximum uncertainty 

■>u 


of a factor of 2 in the calculated value of Q, n /Q n then for E = 0, s 

4-8 x 2.2’ c 

should be smaller than 2.5 0 as we have seen above (Figure 8), s = 4 with 

E = o gives a yery poor fit to the observed pulse-height spectrum, 
c * 

Finally, w'p can estimate the product of the ambient hydrogen 

•\ 

\ 

density and the 'total number of protons of energies greater than 30 MeV, 

nNp(>30 MeV). From the observed 4 to 8 MeV gamma-ray flux of 0.19 
2 1 43 3 

photons cm sec , nN (>30) =* 6xl0 4 cm , for the nominal cases given above 

P 

(s = 1.9, E^ = o and s - 2.5, E c = 30 MeV/amu). However, since the 

E c = 100 MeV/amu case can also give a fit to the data, nN^(>30 MeV) could be 
44 -3 

as large as 2x10 cm . These numbers are in good agreement with those given 
by Ramaty et al. (1975) who have calculated nN^ on the basis of the narrow 
lines alone. 

In summary, the results of sections IV. 1 and IV. 2 are that flat particle 
spectra give a much better fit to the observed gamma-ray spectrum that do 
steep spectra, that steep spectra produce more total radiation in the 4 to 
8 MeV band than observed, and that the entire observed radiation in the 
4 to 8 MeV band is probably of nuclear origin in the form of broad and 
narrow gamma-ray lines . 

3. Upper Limits on the Relativistic Electron-to-Proton Ratio in Flares 
The implication of the last conclusion given above is that an upper limit 
can be set on the ratio of genuine continuum to line (both broad and 
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narrow) emissions in the 4 to 8 MeV energy band. This limit, in turn, sets 
limits on the proton-to-electron ratio. In the discussion below, we take 
into account bremsstrahlung only. The inclusion of other radiation 
mechanisms (which, however, make only very minor contributions) will only 
lower the deduced e/p limit. 

In Figure 13 we show the ratio of the differential proton number to the 
differential electron number at 10 MeV, N^(10) /N e (10) , derived so that 
the total amount of bremsstrahlung in the 4 to 7 MeV range equals Q. . 

It is also assumed that the ambient densities in the interaction regions 
of the nuclei and the electrons are the same. As can be seen, for a wide 
range of spectral parameters, solar flares accelerate protons to energies 
around 10 MeV much more efficiently than they accelerate electrons to 
similar energies,, Bai and Ramaty (1976) have arrived at the same conclusion. 
Furthermore, these authors have argued that relativistic electrons .and nuclei 
of energies 1 MeV are accelerated in solar flares by a different 
mechanism than that which accelerates lower energy (^100 keV) particles. The 
flare mechanism which accelerates particles to high energies must accelerate 
protons more efficiently than electrons. We should mention that the 
accelerator of galactic cosmic rays should have the same property: as for 
solar flare particles, the ratio of protons to electrons in galactic 
cosmic rays is around 100 (e.g. Meyer, Ramaty and Webber 1974). 

V. Conclusions 

We have evaluated in detail the production of nuclear gamma rays in 
the 4 to 8 MeV region,, The radiation in this energy band is due to both 
narrow lines, produced in the interaction of energetic protons and alpha 
particles with ambient material, mostly C and 0, and to broad lines 
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resulting from the interactions of energetic C and 0 nuclei with ambient 
H and He. 

We have studied in detail the nuclear processes which lead to line 

emissions „ We find that the 4 to 8 MeV band has the unique property that 

it contains almost all the nuclear radiation of the abundant elements G 

and 0, but does not receive significant contributions from other nuclei. 

12 11 

This band shows strong lines at 4.44 MeV from C and B, at 6.13 MeV 

from 16 0, at 6.32 MeV from ^N, and at ~5.2 MeV from and ^"*0. The 

6.13 MeV line is produced by direct excitation reactions, the 4.44 MeV 

line is produced by both direct excitation and spallation reactions, and 

the 6.32 and 5.2 MeV lines are due to spallation reactions only. Spallation 

reactions, in general, make a strong contribution to the gamma-ray spectrum 

only if the charged particles have flat energy spectra. 

The FWHM's of the narrow lines are almost independent of the energies 

of the primary particles. The FWHM's of the 4.44-MeV and 6.13-MeV lines 

are about 45 keV and 65 keV respectively, if excited by protons, and 

larger by about a factor of 2, if excited by alpha particles. Since for 

steep particle spectra, most of the gamma rays are produced by alpha 

particle induced reactions, the widths of the lines for such spectra are 
significantly larger than the widths obtained for flat spectra. 

Thus, the widths of the lines as well as the presence or absence of 

lines from spallation reactions could give information on the energy 

spectrum of the accelerated particles. 

The gamma-ray spectrum is also affected by the angular distribution of 

the accelerated charged particles. For example, if the particles have 

velocity vectors confined to a hemisphere either towards or away from 

the observer, the narrow lines are shifted upwards or downwards in energy 0 
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The shift of the mean energy of the 6.13-MeV line, however, is only about 
15 keV, a value too small to be resolved by present detectors. The shift 
of the 4 0 44-MeV line is even smaller. On the other hand, the broad lines 
undergo much larger shifts which lead to modifications in the gamma-ray 
spectra that are observable even with presently available instrumentation. 

We have folded the calculated gamma-ray spectra through the response 
function of the University of New Hampshire gamma-ray detector flown on 0S0-7 . 
By comparing the calculated and observed pulse-height spectra, we conclude 
that the entire 4 to 8 MeV spectrum can be attributed to nuclear radiation, 
that flat particle spectra produce gamma-ray spectra which fit the data much 
better than those produced by steep spectra, and that the data is better 
fit if the particles have isotropic or upward-moving angular distributions 
than if the particles' velocity vectors are confined to the downward 
hemisphere. Since the 1972, August 4 flare was close to the sub-Earth 
point on the Sun, the upward or downward hemispheres correspond to those 
toward or away from the observer. 

Information on the spectrum of the accelerated particles can also be 

obtained by comparing the total nuclear radiation in the 4 to 8 MeV region 

with the 2.2-MeV line intensity observed from solar flares. We find that 

the spectral index of the differential particle number in the 10 to 100 MeV/amu 

range cannot be larger than 2.5, since for steeper spectra the 4 to 8 MeV 

band would contain more radiation than observed. This result takes into 

account uncertainties due to the composition of both the accelerated 

3 

particles and the ambient medium, and the abundance of He in the photo- 
sphere. It is not possible, however, to set a lower limit on the spectral 

index „ Even though very flat spectra produce an overabundance of 
E 

'f\ 
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neutrons relative to 4 to 8 MeV photons, this overabundance does not 
necessarily lead to an overabundance of 2.2-MeV gamma rays, since the 
neutron- to-gamma-ray conversion efficiency decreases with increasing upward 
anisotropy of the accelerated particles. .This anisotropy could be studied 
by measuring simultaneously the fluxes of 2.2-MeV gamma-rays and solar neutrons. 

The fact that essentially all the observed radiation in the 4 to 8 MeV 
region from the solar flare of 1972, August 4 is of nuclear origin, sets 
upper limits on the electron-to-proton ratio in flares. We find that 
at 10 MeV, the ratio between the differential number of electrons 
to protons cannot exceed about 0.05. 

Ac knowl edgmen t s 

We wish to thank Professor E. L. Chupp for providing the 0S0-7 
gamma ray data. Part of the research of A.N.S. was supported by Goddard 
contract NAS5-11276 at the University of New Hampshire. 



-30- 


VIo Appendix 

It is well known that when monoenergetic photons in the 
energy range from about 0.2 to 20 MeV are incident on a NaI(T£) 
detector, the observed pulse-height spectrum consists of a photo- 
peak, Compton continuum and, above 1.022 MeV of two annihalation 
gamma-ray escape peaks. It is often desirable to convert the 
observed Nal pulse-height spectrum into the true gamma-ray spectrum 
or visa versa. The problem of converting a measured pulse-height 
spectrum into an incident photon spectrum requires detailed information on 
the detector response function. 

The response function of the University of New Hampshire (UNH) 
gamma-ray detector was measured for gamma-ray energies of 0.393, 

0.662, 1.12, 1.37 and 2.75 MeV (Higbie et al. 1973). Above 1 MeV the 
photopeak sensitivity of the detector can be approximated by the 
relation 


—0 9 

S , = 10.1 6 * counts cm /photon, (16) 

ph 


where e is the energy of the photon in MeV. We have extrapolated 
this relation to the 4-8 MeV range because there is no experimental 
data above 2.75 MeV. When a flux F of gamma rays of energy e is 
incident on the detector for a time t it causes a photopeak in the 
observed spectrum; the number of events in this peak is given by 


N . = S . F t. 
pn pn 


(17) 


The UNH gamma -ray detector is shielded by a thick CsI(Na) shield. 
The anticoincidence action of the shield reduces the contributions 
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of the Compton continuum and the annihilation gamma-ray escape 
peaks by a factor about 1.5 0 As a result of the suppression of 
the escape peaks, the area under the second escape peak in the 
energy range 4 to 8 MeV is negligible. For the UNH gamma-ray detector 
the ratio between the area under the first escape peak and the 
photopeak was measured at 2.75 and 6.13 MeV. Experimental data 
(Heath 1964) and the Monte-Carlo calculations of Berger and Seltzer 
(1972) show that this ratio is a linear function of the gamma-ray 
energy in the energy range 2 to 8 MeV. We have, therefore, used the 
following relation in the energy range 4 to 8 MeV, 

N /N , = 0.08e - 0ol35, (18) 

es ph * 

where N ' is the number of counts in the first escape peak. 

6S 

Measured pulse-height spectra for gamma-ray energies of 0.66, 

1.12, 1.37 and 2.75 MeV were used to calculate the experimental 

values of the peak-to-total ratio 0 Here is the total 

number of counts in the pulse-height spectrum. In the determination 

of we have approximated the Compton tail and the 'valley fill' 

function by a single rectangular function. Since no experimental 

data is available above 2.75 MeV, we have used the Monte-Carlo 

calculations of Giannini et al. (1970) to construct the peak-to- 

total ratio curve above 3 MeV. Giannini et al (1970) have calculated 

the photopeak and the total detection sensitivities for a 7.6cm x 7.6cm Nal 

crystal for gamma-ray energies up to 10 MeV. The Monte-Carlo calculations of 

Giannini et al. (1970) were chosen because their calculations give peak-to-total 

ratios which are in good agreement with the experimental data (e a g. 

Heath 1964, Misra and Sadasivan 1969) . Above 3 MeV, the peak-to 
total ratio for the UNH gamma-ray detector can be represented by 
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the relation. 


V /n t " °- 69 e 


-0.80 


(19) 


Knowing N N e s* * nc * ^T* t ^ ie nutn ^ er counts in the Compton 
continuum is given by 

N = N - N , - N . (20) 

c T ph es 

The relationship between the energy of a gamma ray and the 

channel number n of the pulse height analyser is given by a 

quadratic formula (P. P. Dunphy, private communication 1973) 

8 = K(n+n Q ) 2 , (21) 

where from ground calibration data n =80.2. The value of 

o 

K can vary due to gain changes in the detector. For the 1972, 

—4 

August 4 flare period the mean value of K was found to 0.393x10 

The procedure adopted to convert the calculated photon spectrum 
in the energy range 4-8 MeV to the pulse-height spectrum is as 
follows. The energy range 4-8 MeV was divided into bins of 200 keV 
width each 0 This bin width is chosen to retain the instrumental 
resolution width of the detector which is about 140 keV at 5 MeV. 

All the photons falling in a given energy bin are considered to have 
the mid-bin energy of chat bin. Starting from the highest energy 
bin, the contribution of the Compton continuum associated with the 
counts in this bin was determined for all the bins lower than 
this bin. This procedure was repeated for the second bin and so on 
down to the last bin. The total Compton contribution for each bin 
was then determined by summing up the individual contribution from 
bins higher than that bin. Similarly starting from the highest 
energy bin and working down, the area under the first escape peak 
was determined and added to the bin shifted down in energy by 
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0.51 MeV from the contributing bin. Finally, the expected number 
of counts in each bin was determined by adding the photopeak, Compton 
and first escape peak contributions. 

It is clear from the foregoing discussion that the counting 
rate, due to a flux of photons from a gamma-ray spectrum incident 
on the detector, consists of the sum of contributions from photopeak, 
first escape peaak and Compton continuum. Using the folding procedure 
described above, we have converted several calculated gamma-ray spectra 
into pulse-height spectra. From the comparison of these spectra 
we find that the relation between the incident flux of 4 to 8 MeV 
photons (F^ g) and the counting rate (C^) can be expressed as, 

C R " R F F 4-8 " °* 106 F 4-8- (22 > 

This relation depends only weakly on the spectral shape of the 

accelerated charged particles. For the particle spectra discussed 

in this paper, the value of C does not change by more than 5% 

of the value given in equation (22) . 

The observed pulse-height spectrum for the 1972, August 4 event 

was accumulated during the time interval of 91.6 sec. The total 

number of counts observed in this time interval was 82+9. Since 

2 

the area of the detector is 45cm , F, 0 = 0.19+0.02. This is the 

* 4-8 — 

number that we used in Section IV. 2. 
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Table 1 


Prompt gamma ray lines in the energy range 4 to 8 MeV 


Photon Energy 
(MeV) 

4.438 


4.443 

5.180 

5.241 

5.270 

5.298 

6.129 

6.176 


6.322 

6.337 

6.478 

6.741 

6.791 

6.917 


Emission 

Mechanism 

Production 

Processes 

i2 C* 4 ,439 -*g .s . 

12 C(p,p' ) i3 C* 
13 C(a,c/) 13 C* 
1B 0(p, ) 13 C* 

i4 n(p; > i2 c* 

lig*4 *444 

13 C(p,2p) ll B* 
l3 C( a , ) 11 B* 

lB 0(p, 

l4 NCp ) i:l B* 

lo 0 *6-lal_ gsS . 

i £ 0 * 6,34s - g .s. 

16 0(p,pn) 1B 0* 

16 0(a, ) l5 o* 

lb N *S.y71_^ gs> 

* 5 N* 5,299 -g. S . 

i6 0(p,2p) 15 N* 
16 0 (a, ) 15 N* 

16q*6 . 

16 0(p,p / ) 1B 0* 

lb 0(a,or / ) ls 0* 

i 5 o* 6 * 177 -»g .s . 

16 0(p,pn) ib 0* 
l6 0(a ) 15 0* 

lE N* 6 * 324 -*g.s. 

i6 0(p,2p) lB N* 

l6 0( a , ) 1& N* 

ll c * 6 * 339 -» g . s . 

12 C(p,pn) 11 C* 

ia C *6 *480_ g>s> 

13 C(p,pn) 11 C* 

ll B *6.743^ g<s> 

12 C(p,2p) 11 B* 

iX B *B-793^g eS> 

’• 2 C(p,2P) 11 B* 

16 0 *6 'String >g# 

iB 0(p,p' ) lB 0* 

ie 0(a,a' ) i6 0* 

16 0 *7 -liS^g.s. 

l6 0(p,p' ) 16 o* 
16 0 a, a' ) ie 0* 


7.117 
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Table 2 

Line Widths of 4.44 and 6.13 MeV Lines 


S E c FWHM (keV) 

(MeV/amu) 4.44 MeV 6.13 MeV 

4 0 110 140 

20 45 65 

2 30 45 65 

2 100 50 120 
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Table 3 


-1 -1 

Selected 4-0 MeV Production Rates (photons sec Hatom ) and 
Neutron Production Rates (neutrons sec ^ Hatom 



He/H = 

0.07 

He/H - 

0.1 


or/p 

Q 4-8 

^n 

Q 4-8 



0.02 

2.37xlO~ 18 

1.18xl0" 16 

2.38x10 AO 

1.39xl0 -16 


0.05 

3.44x10" 18 

1.43xl0 -16 

3.47xl0“ 18 

1.63xl0 -16 ' 

S, S - 1.5 
{ E = 30 MeV/amu 

0.1 

5.23xl0" 18 

1. 83xl0 -16 

5. 29xl0 -18 

2. 03xl0 -18 - 

s c 

0.02 

1.66xl0" 15 

0.92xl0 -15 

-15 

1.97x10 

-15 ' 

1.14x10 


0.05 

0.1 

3.68xl0~ 15 
7 .04xl0~ 15 

2 . 12xl0 -15 
4 . llxlO -15 

4.55xl0 -15 

8.57xlO -15 

2.65xl0 -15 1 
-15 

5.16x10 , 

> S = 5 

E = 0 
c 
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Figure Captions 

12 

1. Production cross sections for the 4.44 MeV line from C. Closed 
data points - proton induced reactions; open point - alpha 
particle induced reactions. 

2. Production cross sections for gamma-ray lines from spallation 

16 

reactions on 0. Closed points - 4.44 MeV gamma rays; open 
points - 5.2 MeV gamma rays. The dashed-dotted line is the sum 
of the cross sections for the production of the 6.322 and 6.176 MeV 
lines. 

16 

3. Cross sections for line production in the 6. 1-6. 3 MeV range from 0. 

4. Production cross sections for the 2.74, 6.92 and 7.12 MeV gamma-ray 
lines from ^0 o 

5. Production rates of narrow gamma-ray lines in the 4 to 8 MeV region. 

The sum of the production rates of the narrow lines due to proton 
(q^) and alpha (q^) particle interactions are also shown in this 
figure. 

6. Calculated gamma-ray spectra resulting from accelerated charged 
particles with isotropic angular distribution. 

7. Calculated gamma-ray spectra resulting from anisotropic particle 

distributions with energy spectrum characterized by s = 2 

and E =30 MeV/amu. 
c 

8. Comparison of the observed pulse-height spectrum from the flare of 
1972, August 4 with the calculated pulse-height spectra obtained 

from the gamma-ray spectra shown in Figure 6. We limit the comparison 
to the 4-7 MeV range because there were no observed counts above 7 MeV. 
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9o Comparison of the observed pulse-height spectrum from the flare of 
1972, August 4, with the calculated pulse-height spectra resulting 
from the gamma-ray spectra shown in Figure 7. 

10. g as a function of s and E c , where g is the production rate 
of photons in the 4 to 8 MeV energy interval. These production 
rates are normalized so that for all values of E^, including 

E £ = 100 MeV/amu, the number of protons with energies greater than 
30 MeV is unity. 

11. Q as a function of s and E , where Q is the production rate of 

n c n 

of neutrons. The normalization of the rates is the same as that 
in Figure 10. 

12. The Q. 0 /Q 0 0 ratio, where Q 0 0 is the production rate of 2.2-MeV 
gamma rays from neutron capture. The shaded area represent the 
observed ratio of photons in the 4 to 8 MeV range to the observed 
photons in the 2.2 MeV line. 

13. Np(10) /N^ClO) as a function of the spectral index (Se) of the 

accelerated electrons. N (10) and N (10) are the differential 

P e 

numbers of protons and electrons at 10 MeV respectively. 


$ Zobel et at (1968) 
$ Alard et al (1974) 
$ Zobel et al (1968) 


X. \' 


\ \\ x 

V W 


l2 C(p,p')' 2 C* 44 Y s n v 

l2 C(a,a')' 2 C* 4 ' 4 ^ 

l2 C(a, a') l2 C* 4 - 4 + ,2 C(a,a P ) "b* 4 - 4 
l2 C(p,p') l2 C* 4 - 4 + l2 C(p,2p)"B* 4 - 4 













10 MeV/amu 



H : C : N : 0 = 10*: 3.7: 1.2: 6.8 

N a (E) = 0.1 N P (E) 

m t c \ /const E- E c 
Np(E)oc | E . S E> ^ 

N P (>30 MeV) = I 
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